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Abstract 



We measure the fractal structure of four dimensional simplicial quantum gravity 
by identifying so-called baby universes. This allows an easy determination of the 
critical exponent 7 connected to the entropy of four- dimensional manifolds. 
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1 Introduction 



A simple way to characterize the fractal structure of 2d quantum gravity was recently 
introduced in ref. Using very few assumptions it was possible to prove that the 
functional behaviour of the number of two-dimensional surfaces with area N 2 and 
coupled to matter with central charge c < 1 

AT{N 2 ) ~ iV 2 7(c) " 3 e^ N2 (1) 

is directly related to the fractal structure of 2d quantum gravity. In particular one 
can determine 7(c) by simply counting the number of so-called "minimum bottleneck 
universes", abbreviated "minbu's", living on a typical surface in the ensemble of 
surfaces of fixed total area. This number is given by 

M Ni {V 2 ) ~ N 2 V^ C) ~ 2 , (2) 

where N 2 is the total area of the surface and V 2 <C N 2 is the area of the minbu's 
being counted. A "baby universe" of area V 2 is any simply connected region of 
the surface of area V 2 and boundary length I such that V 2 ^> I 2 , and a minbu is a 
baby universe whose / is the smallest possible length consistent with the ultraviolet 
cutoff. This boundary, along which the minbu is connected to the remainder of 
the (parent) surface is called a minimal bottleneck. The only assumption going 
into the derivation of (0) is the distribution ([[]) and the approach is well suited 
for dynamically triangulated surfaces 0, |3|, f|, |5|. Here the area of a region can be 
identified with the number of triangles in the region and a minimal bottleneck is 
a simple closed path consisting of three links which partitions the surface into two 
parts - the minbu containing V 2 ^> 1 triangles, and the parent containing N 2 — V 2 
triangles. 

Until now the measurement of 7 by numerical simulations has been a rather 
painful procedure since one had to perform a whole sequence of simulations (or a 
so-called grand canonical simulation where the number of triangles changes, which 
is also an unpleasant task) for different areas N 2 and in the end fit the measured 
quantities to a formula related to ([!]). Following the suggestion in |] one can now 
extract 7 from a single simulation using ([|). This is neater also because one does 
not have to subtract out the leading e^° N2 piece while making the fit; (^]) contains 
only the universal power law piece. This method also works in practise. Indeed, in 
a recent numerical simulation it was possible to extract 7 for pure 2d gravity with 
high precision (less than 1%) 0. 
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2 The model 



A model of 4d quantum gravity which generalizes the discretized 2d approach was 
suggested a year ago []7], ||. The partition function is given by 

Z(k 2 ,k 4 )= Y,e- K ^ + ^ (3) 
TeT 

where the sum is over triangulations T in a suitable class of triangulations T. The 
quantity N4 denotes the number of 4-simplexes in the triangulation and N 2 the 
number of triangles. The coupling constant k 2 is inversely proportional to the bare 
gravitational coupling constant, while k 4 is related to the bare cosmological constant. 
The most important restriction to be imposed on T is that of a fixed topology. If 
we allow an unrestricted summation over all topologies in @ the partition function 
is divergent ||. In the following we will always restrict ourselves to considering only 
manifolds with the topology of S 4 . 

Z(k 2 , K4) is the grand canonical partition function. It is defined in a region 
K4 > kI(k 2 ) i n the (k 2 , ^4) coupling constant plane. The only way in which we can 
hope to obtain a continuum limit is by letting k 4 approach ^4(^2) from above. This 
tentative continuum limit depends only on one coupling constant k 2 . We can write 
© as 

Z(k 2 ,k,) =Y,Z(K 2 ,N 4 )e- K * N \ (4) 

N4, 

Z(k 2 , N4) is the canonical partition function where iV 4 is kept fixed. Then we have 
in practise only one coupling constant, k 2 , and the aspects of gravity which do not 
involve the fluctuation of the total volume of the universe can be addressed in the 
limit of large iV 4 . However it is actually possible to extract the critical exponent 
for volume fluctuations from Z{k 2 ,N 4 ). Let us assume that the canonical partition 
function has the form: 

Z(k 2 , N a ) = N2 (K2) ~ 3 e^ R * )Ni ■ (1 + 0(1/N 4 )) (5) 

When K4 is close to k1(k 2 ) we can approximate (^) with 

Z( K2 , .4) ~ analytic + ^ _ ^ ))l(K2) „ 2 (6) 

where "analytic" means possible analytic terms of the form (k 4 — n c A ) n 1 < n < 
2 — 7(^2). The entropy exponent of the number of 4d simplicial manifolds, 7, 
determines the volume fluctuations since we have 

AT 2\ , AT \2 d 2 \nZ(K 2 ,K 4 ) C(k 2 ) 

N * > ~ {N4) = M ana ytlc + (.4 - kim)^) (7) 
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if we assume 7(^2) < 2. 

We will try to extract the critical exponent 7(^2) by numerical simulations. At 
this point it is important to note that it is by no means obvious that the exponent 7 
exists. Equation (|5p is an ansatz, inspired from 2d quantum gravity where we know 
that a similar formula holds when the central charge of matter c < 1. In the case 
of c = 1 there are logarithmic corrections, and the form is not known for c > 1. In 
3d simplicial gravity it was shown [15], [17], [l8j that a more likely form of Z(k±, N 3 ) 



(the 3d analogue of @) is 

C{k x ) 



Z(ki,N 3 ) ~ exp 



k%( Ki )N 3 1 - 



K 
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and a possible power law correction would be subleading. However, 3d simplicial 
gravity seems to differ from 4d simplicial gravity in many respects (0, [8], |16|), and it 
was clear from the fine tuning process K4 — > kI(k2) in the Monte Carlo simulations 
that the corrections to the leading term exp^iV^) were less severe in 4d than in 3d. 
It is therefore indeed possible that the finite size corrections in 4d are power-like, as 
in 2d simplicial gravity, rather than of the exponential form (pi) encountered in 3d. 

Let us now define the minimal bottleneck baby universes, which we, following 
will denote "minbu's". On a 2d triangulated manifold one can check whether we 
have a minbu in the following simple way: pick a link and two neighbouring links, 
attached to the two vertices of the first link. If the two links have a vertex in common 
the three links form a closed loop. In general this loop will be trivial in the sense 
that it will just be the boundary of one of the two triangles in the surface which 
contains our starting link. However, it might turn out that there is no triangle in the 
surface which has the closed loop as its boundary. If we cut the surface along such 
a closed loop we will have separated it into two disconnected open surfaces each one 
having the closed loop as its boundary, provided the topology of the surface is that 
of a sphere (which we will assume). Both components will have spherical topology 
when closed by adding the triangle which has the loop as boundary. The smallest 
of the spheres will be a minbu, the largest the "mother universe" . 

This construction can immediately be generalized to 4d triangulated manifolds. 
We check for minbu's as follows: pick a 3-simplex (a tetrahedron) and in this the 
four 2-simplexes (triangles) which constitute its boundary. Identify for each triangle 
all the 3-simplexes in the manifold which have this triangle as their boundary. We 
now have four groups of 3-simplexes and all 3-simplexes have one vertex which does 
not belong to the triangle from which the 3-simplex was constructed. Pick now a 
3-simplex from each of the 4 groups and check whether their 4 free vertices coincide. 
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If that is the case the four 3-simplexes together with the original 3-simplex can be 
thought of as the (closed) boundary of a 4-simplex. In general the 4-simplex found 
by this procedure will just be one of the two 4-simplexes which shared the original 

3- simplex. However, if that is not the case we will, if the topology of our original 

4- manifold is that of a sphere (which we assume as usual), separate it into two 
non-trivial components if we cut along the closed piecewise linear 3d manifold built 
out of the five 3-simplexes. Both components have this closed 3-manifold as their 
boundary and both components will have the topology of S* 4 if we close them by 
adding to each of them a 4-simplex in such a way that their boundary is identified 
with the boundary of the 4-simplex, just with opposite orientation. The smallest of 
the 4-spheres will be called the minbu, the largest the "mother" , in agreement with 
the 2d notation. 

It is almost clear that the above description can be turned into an efficient nu- 
merical algorithm (after, admittedly, some pain with double counting etc.) provided 
one is given the coincidence matrix of the triangulated 4-manifold. 

Let us now argue that the counting arguments in extend to the 4d triangulated 
manifolds as well. Let us by Z'(k2, N4) denote the canonical partition function for 
4-manifolds consisting of N4 4-simplices where one 4-simplex is marked. Generically 
we will get a different manifold for each mark so for large N4, where accidental 
symmetries are expected to play no role, we get Z\kz,N^) ~ N±Z(k2, N4). Such a 
marked manifold can also be viewed as a manifold with (iV 4 — 1) 4-simplexes and 
a minimal boundary of the kind associated with minbu's (by removing the interior 
of the marked 4-simplex). A moments reflection will convince the reader that the 
average number of minbu's of volume V4 on 4-manifolds of total volume N4 will be 
given by 

60 

where 60 is the number of ways one can glue the two boundaries of the minbu and 
the mother together with opposite orientation of the boundary. 

If we assume the canonical partition function is given by (||), we get 

■Afatei, v A ) ~ c{k 2 ) n a v2 {K2) - 2 (10) 



The scene is now set for a numerical determination of the number of minbu's since it 
is already by now |7], ||, [H], [TTJ, [L^, [13], [L4| standard how to generate by Monte Carlo 
simulations the class of triangulated 4-manifolds corresponding to the partition func- 
tion (0). It is also well known || how one can effectively stay in the neighbourhood 
of a given volume, iV 4 , even if one is using the grand canonical partition function 
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(|3p, as one is forced to by ergodicity requirements in 4d, contrary to what is the 
case in 2d. 



3 Numerical results 

By Monte Carlo simulations we have generated a number of 4d-manifolds according 
to the distribution dictated by the partition function (H). For details about this 
procedure we refer e.g. to [^, |T2||. We have considered k 2 = 0.0, 0.8, 1.0, 1.1 and 
1.4 and N A = 4000, 9000, 160000. For each value of k 2 and for each value of N 4 we 
have generated approximately 1000 configurations. 

Let us at this point review some of the results of For k 2 ~ 1.0 — 1.1 we 
observed a transition in geometry from a highly connected phase with a seemingly 
large Hausdorff dimension to a phase with an elongated, almost one-dimensional 
geometry. This is why our main efforts have been concerned with this region of 
coupling constant space. 

The thermalization time is short in the highly connected phase, but quite long 
for k 2 > 1.0. In the case of iV 4 = 16000 we have for these values of the coupling 
constant used 10000 sweeps^ for thermalization and performed measurements after 
each successive fifty sweeps (after each successive hundred sweeps for k 2 = 1.4). The 
result of the measurements of minbu's is shown in fig. 1 for N4 = 16000. We have 



plotted the number of minbu's on a log- log scale since the formula (10) suggests the 
dependence: 

logA6v 4 (K 2 , V 4 ) = Const(«; 2 ) + ( 7 (k 2 ) - 2) \ogV A (11) 

The straight lines are the result of a least x 2 fit- In fig. 2 we have shown the results 
for fc 2 = 1.1 and different values of N± and finally in fig. 3 we have shown 7(^2) as 
extracted from fig. 1. 

A few comments should be made. We have only taken into account minbu's with a 
volume V4 > 9. For small values of V4 the numbers A/jv 4 (^2, V&) fall into two classes 
according to whether V4 = 4n — 1 or V4 = 4n + 1. This is a clear finite size effect. 
We have only included in our analysis the last class of numbers as they seem to fit 
(|TT|) all the way down to V4 = 9. For large volumes V4 of the minbu's the results are 
average values for volumes in the neighbourhood of the plotted V4 since the statistics 
for a specific large volume is not good. The errorbars resulting from the binning 
procedure are smaller than the symbols used in the figures. 



2 We have also made a few runs for N4 — 32000 and k-i = 1.0. 
3 By a sweep we mean N± accepted updatings. 
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In general we see a clear deviation from (O) when MiV^) < 0.1 and when 
«2 < in this case the number of minbu's begins to drop faster than given by 
(11). We have not shown these data, since they have bad statistics, and we do not 
know whether the above mentioned behaviour reflects that ([□]) is not really valid for 
large V4 or rather that the statistics is just not good enough (there is typically less 
than one such large minbu per universe) for the volumes of manifolds we consider. 
The message we get by comparing different AVs does not allow us at the moment to 
answer this question. For k 2 larger than 1.1 the situation is different: The number of 
baby universes is much larger and even if the number of very large baby universes is 
small and errorbars large their numbers seem not to drop below numbers predicted 

by O- 

It is seen that the curve for k 2 — 1-4 bends at V4 ~ 30 — 40 to a different slope. 
We have used the part with the larger values of V4 to extract 7(^2) ■ Such a bending 
is not observed for k 2 < 1.0, i.e. in the phase with a highly connected geometry. A 
close look at the data for k 2 = 1.1 reveals a slight tendency to such a bending, and 
since k 2 = 1.1 is just at the borderline of the transition to the elongated geometry 
the bending is probably related to the fact that the very fractal structure observed 
in this phase is only unfolded for large volumes V4. 



4 Discussion 



With the above reservation concerning the interpretation of the numerical data, we 
still consider the results as somewhat remarkable. It seems possible by numerical 
simulations to study in detail the fractal structure of quantum gravity. Indeed, our 
data contains a lot of information about the branching of the geometry which will 



be published in a more extensive report elsewhere ||19| . Let us just here mention 
that the two phases of 4d simplicial quantum gravity observed in former works ap- 
pear quite clearly in the minbu measurements. In the phase with highly connected 
geometry there is one "big mother" which contains the major fraction of the volume 
(from 95% to 80% depending on « 2 )- In the other phase of very elongated geometry 
this is not so and there seems genuine democracy in "mother size" . Nevertheless the 
7(^2) extracted seems perfectly smooth when we pass the transition in geometry, as 
is apparent from fig. 3. It is interesting to note that a similar smoothness has been 



observed in the numerical simulations of 2d quantum gravity [^, |20| when we pass 
the c = 1 barrier. It points to the fact that subleading corrections to (|TJ) might play 
a very important geometrical role, and it might be possible to understand analyti- 
cally the interplay between the subleading terms in ([I]) and the fractal geometry of 
quantum gravity. The seed to such an analysis is already to be found in M. 
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We would like to stress that the results here, especially for k% < 1.1, are obtained 
for quite small baby universes and we see some deviations for larger volumes, as 
mentioned above. It would clearly be desirable to let the babies grow, but as is 
well know from biology this requires time, in our case computer time. Work in this 
direction is in progress. It is nevertheless somewhat intriguing that the transition 
between the two types of geometry takes place at a value of 7 ~ —0.5 — 0.0, i.e. 
the same values of 7 which are of main interest in 2d gravity. Note further that 
the results for 7 seem internally consistent in the sense the value in the phase with 
almost linear geometry (Hausdorff dimension not larger that two) seems close to the 
value 7 = 1/2 for the so-called branched polymers. In addition, assuming that the 
extraction of 7 for K2 < 1.1 will survive the test of large babies, one could wonder if 
this reflects that the situation is like in 2d gravity where we have a 7(c) depending 
on the central charge. In our case k 2 would then play a role somewhat similar to the 
central charge. The range of 7's is in agreement with such a picture. An interesting 
continuum approach somewhat related to such an interpretation can be found in a 
recent paper by Antoniadis, Mazur and Mottola [21]], where the authors develop a 
quantum theory for the conformal mode of gravity. 
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Figure Captions 



Fig. 1 The number of minbu's for iV 4 = 16000 and k 2 = 0.0(D), 0.8(0), 1-0(A), 
1.1 (+) and 1.4(x). The straight lines represent the best fits to ([10|). 

Fig. 2 The number of mimbues for k 2 = 1.0 and iV 4 = 4000(D), 9000(O), 16000(A) 
and 32000 (+), normalized such that Mn a {V) is multiplied by 16000/A^. (log- 
log scale) 

Fig. 3 7(^2) as extracted from fig.l. 
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